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■ Abstract. We present wide-field spectroscopy of globular clusters around the Leo I group galaxies NGC 3379 and NGC 3384 

using the FLAMES multi-fibre instrument at the VLT. We obtain accurate radial velocities for 42 globular clusters (GCs) in total, 
30 for GCs around the elliptical NGC 3379, eight around the lenticular NGC 3384, and four which may be associated with either 
. galaxy. These data are notable for their large radial range extending from 0^7 to 14^5 (2 to 42 kpc) from the centre of NGC 3379, 

' and small velocity uncertainties of about 10 km s"' . We combine our sample of 30 radial velocities for globular clusters around 

, NGC 3379 with 8 additional GC velocities from the literature, and find a projected velocity dispersion of CTp = 175^2? kms"' 

alR < 5' and <Tp = 147+3^ alR > 5'. These velocity dispersions are consistent with a dark matter halo around NGC 3379 with a 
concentration in the range expected from a ACDM cosmological model and a total mass of a; 6 x 10"Mo. Such a model is also 
consistent with the stellar velocity dispersion at small radii and the rotation of the Hi ring at large radii, and has a (M/L)b that 
^ increases by a factor of five from several kpc to 100 kpc. Our velocity dispersion for the globular cluster system of NGC 3379 

, is somewhat higher than that found for the planetary nebulae (PNe) in the inner region covered by the PN data, and we discuss 

possible reasons for this difference. For NGC 3384, we find the GC system has a rotation signature broadly similar to that seen 
in other kinematic probes of this SBO galaxy. This suggests that significant rotation may not be unusual in the GC systems of 
, disc galaxies. 

^. 
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1. Introduction 

The current paradigm for galaxy formation is that galaxies form 
from baryons that cool within dark matter halos. This paradigm 
is supported by the abundant observational evidence for dark 
matter in the Universe on many scales, ranging from clusters 
of galaxies to the halos of individual spiral and dwarf galaxies. 
Ideally, one would like to test for the presence of dark matter 
halos around early-type galaxies as well, and to study the prop- 
erties of these halos. 

Determining the halo masses and mass profiles of individ- 
ual elliptical galaxies has proven to be challenging, primar- 
ily because they lack a readily observed dynamical tracer at 
large radii like the H i gas of spiral galaxies. Dedicated spectro- 
scopic studies of the integrated light have been limited to about 

Send offprint requests to: S. Zepf, e-mail: zepf @pa . msu . edu 

* Based on observations collected at the very large telescope (VLT), 
European Southern Observatory, Chile, under programme 71.B-0537. 

** Tables IXTI IA.2l and lA.3l are available in electronic form at CDS 
(ftp cdsarc . u-strasbg . f r) or http : //www . edpsciences . org. 



R < 2Rs (e.g., Kronawitter et al. I2000> . For the most lumi- 
nous and massive ellipticals. X-ray observations of the hot gas 
around these galaxies provide evidence for dark matter (e.g. 
Loewenstein & White 0^55). Studies of gravitational lensing 
also indicate the presence of dark matter halos around ellip- 
tical galaxies. Specifically, individual cases of strong lensing 
(e.g., Treu & Koopmans I2004I I show evidence for dark mat- 
ter at modest radii, while analyses of weak lensing (Wilson 
et aL lMJTt McKay et al. UMU Hoekstra et al lMHt provide 
strong statistical evidence for dark matter halos at large radii 
around early-type galaxies. Similar results for dark matter at 
large-scales around elliptical galaxies are provided by studies 
of satellite galaxies in large surveys (Prada et al. 120031 Brainerd 
I2004I I. However, for ordinary elliptical galaxies, a substantial 
range of profiles from inner to outer regions is allowed, and 
the dependence of halo properties on galaxy luminosity is not 
strongly constrained. Therefore, determining the masses and 
mass profiles of ordinary individual elliptical galaxies is of sig- 
nificant interest. 
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The radial velocities of globular clusters (GCs) and plane- 
tary nebulae around elliptical galaxies provide a powerful way 
to probe the dynamics of their host galaxy. The extended spa- 
tial distribution of GC systems makes them particularly useful 
probes of any dark matter halos around elliptical galaxies. The 
value of these probes has been demonstrated by the clear evi- 
dence for dark matter halos found by studies of the radial ve- 
locities of GCs around the elliptical galaxies NGC 4472 (e.g. 
Zepf et al. UHM Cote et al. EUU^, NG C 4486 (e.g., Cohen 
l20nol Romanowsky & Kochanek l200n, NGC 1399 (Richtler 
et al. l2004t . and NGC 5128 (Peng et al. l2004l . However, there 
are only four galaxies with such high quality data, and with the 
exception of NGC 5128, all are considerably more luminous 
than a typical elliptical galaxy. 

Spectroscopic studies of the planetary nebulae (PNe) can 
also probe the dynamics of elliptical galaxies. Romanowsky et 
al. (2003j, hereafter R03, specifically carried out a study of the 
radial velocities of PNe around three elliptical galaxies of mod- 
erate luminosity. They found that the galaxy with the most data, 
NGC 3379, required little dark matter to the radial limit of their 
sample at 9 kpc (3')- The other two ellipticals, NGC 821 and 
NGC 4494 also appeared similar, but with weaker constraints. 
NGC 3379 (s M105) is in many ways the archetypal ordinary 
elliptical galaxy, as it is the nearest "normal" giant elliptical 
galaxy, with a distance of 10 Mpc (Jensen et al. 120031 see also 
Gregg et al. Bt)04? in the Leo 1 group, and an intermediate lu- 
minosity of Lb = 1.42 xIO'^Lq (Mb = -19.9). Therefore, the 
possibility that its dark matter halo might not be as expected is 
critical to follow up. 

In this paper, we present the results of a study of the dy- 
namics of the GC system around NGC 3379. Studying the 
GC system has two important advantages for understanding 
the halo dynamics of NGC 3379 beyond the obvious utility of 
providing an independent test of the PNe results. Firstly, there 
has long been discussion in the literature about the possibility 
that NGC 3379 has a significantly flattened component viewed 
face-on (Capaccioli et al. ll99H . Subsequent papers have come 
to mixed conclusions about the presence of such a compo- 
nent (e.g., Statler [T994l Statler & Smecker-Hane [T999l Statler 
1200 U . but it is obviously an important consideration for inter- 
preting the radial velocities of objects observed in NGC 3379. 
Here, GCs provide a key advantage, because they are less likely 
than other tracer populations such as PNe to have significant 
contribution from a strong disc component. The second advan- 
tage of GCs is that candidate GCs can be identified in wide- 
field imaging of nearby ellipticals out to very large radii (e.g. 
Rhode & Zepf l200ni2n04l - hereafter RZ04; Harris et al.I^OOTi. 
This large area over which GCs are observed allows the kine- 
matics of the halo to be probed over a very large radial range. 

The plan of the paper is that the spectroscopic observations 
and the determination of the radial velocities are presented in 
Sect. 121 The properties of the GC systems of the Leo 1 group 
galaxies NGC 3379 and NGC 3384 are discussed in Sect.|3l In 
Sect. 13 we present the dynamical analysis of the NGC 3379 
GC system and the implications for the dark matter halo of this 
galaxy. The conclusions are given in Sect.|5l 



2. Observations and data analysis 

2.1. Identification of globular cluster candidates and 
MEDUSA fibre allocation 

The advent of large format images with mosaiced arrays of 
CCDs on 4-m class telescopes allows for the identification of 
candidate GCs over a wide field around galaxies in the lo- 
cal universe. In particular, RZ04 carried out a multi-colour 
imaging study of the GC systems of nearby elliptical galax- 
ies utilising the MOSAIC camera on the Mayall 4 m telescope. 
They identified GC candidates as unresolved objects with BVR 
colours consistent with those of Galactic GCs, taking into ac- 
count the photometric uncertainties (see RZ04 for more de- 
tails). In order to obtain sufficient S IN spectra in a ^ 4-hour 
exposure, we selected candidates brighter than V < 22. To help 
minimise contamination from non-GCs, we also placed a bright 
limit on the sample ofV> 19, which affects only the extremely 
high luminosity tail of the GC population at My ^ -11. Over 
the 37' X 37' MOSAIC field centered on NGC 3379, 196 can- 
didate GCs met these criteria for our spectroscopic follow-up. 

We used FLAMES, the multi-object, wide-field fibre spec- 
trograph of VLT/UT2 (see Pasquini et al. I2004> to efficiently 
obtain spectra of a significant sample of these GCs. The 
GIRAFFE/MEDUSA mode of this instrument offers 130 fi- 
bres over a 25' diameter field of view, making it an excel- 
lent match to the wide-field MOSAIC images (see Fig.^ and 
the extended spatial distribution of the globular cluster sys- 
tems of nearby elliptical galaxies. We ran the FPOSS soft- 
ware to allocate the 130 available fibres to 108 GC candidates, 
twelve sky positions, and ten V ~ 1% stars to be used as tem- 
plates for radial velocity determinations, in a field centered at 
{a,6)j2ooo.o = (10h48'"0P, 12°33'30"). 

2.2. Observations and radial velocity determination 

The FLAMES observations were carried out in service mode 
on VLT/UT2 (Kueyen), with the GIRAFFE/MEDUSA config- 
uration. We chose to use the grating LR4 which has a resolv- 
ing power of A/A(A) ~ 6000, and a central wavelength of 
Ac = 543 1 A that covers many absorption features over 800 A 
(including strong Mg and Fe lines), and sufficient efficiency to 
provide useful spectra (S /N > 5 per pixel) of $ 22 unre- 
solved sources in 3-4 hours. Our observing programme con- 
sisted of five 2595-s exposures taken on separate nights, in 
May-June 2003, under good seeing (~ 0'.'7) conditions. 

After standard reductions to extract the spectra with the 
IRAF HYDRA package, we determined the radial velocities v of 
the objects by cross-correlating their spectrum with stellar tem- 
plates. The cross-correlation function (CCF) peak was fitted by 
the IRAF/rv task f xcor to estimate the velocity relative to the 
template, with a typical accuracy of 5-15 kms"' (each pixel 
represents 0.2 A or 11 km s"' with the LR4 grating). 

The cross-correlation was first done using the ten y ~ 1 8 red 
stars to which we allocated fibres so as to provide stellar tem- 
plates for radial velocity determination. To supplement these 
"simultaneous" templates, we used 100 F5-K2 high signal-to- 
noise spectra - both (sub)dwarfs and giants - extracted from 
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Fig. 1. V^band MOSAIC image of the Leo I group. The image 
covers 37' x 37' and is plotted with north up and east to the 
left. The large circle shows the FLAMES 25' diameter field of 
view, corresponding to 75 kpc at the distance of NGC 3379. 
The objects are noted with various symbols according to their 
FLAMES radial velocities v as follows: squares represent GCs 
(500 < V < 1300 kms"') with larger squares indicating a 
more certain radial velocity, crosses represent objects identi- 
fied as stars {v < 350 kms"') with the size of the symbol cor- 
responding to the confidence in the velocity as above, and tri- 
angles show the 21 objects without a clear cross-correlation 
peak including objects with insufficient S IN and background 
emission-hne galaxies. 



the ELODIE archive' (Moultaka et al. l2004t . The LL_ELODIE 
stellar library is well-suited for GIRAFFE, with a very similar 
dispersion (the scale is also of 0.2 A/pixel, identical to the LR4 
grating) which implies minimal rebinning of both spectra. 

The CCF peak Gaussian fits to all of the templates were 
checked both by determining the 'R^^^-y?A\it of the fit as de- 



fined by Tonry & Davis (I1979> . and additionally by eye. These 
agreed well, and the velocities obtained for a dozen of the 
brightest targets were compared between individual exposure 
spectra, and also showed good agreement to within the errors 
estimated by fxcor. We then adopted the following classifica- 
tion of the quality of FLAMES radial velocity determinations; 

- class A: more than 98% of ELODIE and more than five (out 
of ten) simultaneous templates agree within the errors 6v; 

- class B: more than 80% of ELODIE and > 3 simultaneous 
templates agree within the estimated errors; 

- class C (not used): less than 60% of ELODIE and only 0-2 
simultaneous templates agree. 
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Fig. 2. A histogram of the FLAMES radial velocities for all 
class A (34 GCs and 32 stars) and class B (8 GCs and 4 stars) 
objects. This histogram shows a clear ~ 200 kms"' gap be- 
tween the highest velocity stars and the lowest velocity clusters 
(associated with NGC 3384). 



For one object, gc771, we clearly detected [Ora]/I5007 at a 
radial velocity of 766 kms"'. This both provides interesting 
evidence for a planetary nebula in an extragalactic GC, suggest- 
ing these may not be exceedingly rare (see Minniti & Rejkuba 
ED02i for another example), and also confirmation of the best 
estimate for the radial velocity from the cross-correlation of 
the absorption lines. For all subsequent work, we consider the 
objects with class A and B radial velocities, and defer consid- 
eration of the class C objects until further data are obtained. 

The reliability and accuracy of our radial velocities are 
strongly supported by a comparison with the data of Puzia et 
al. 120041 (hereafter P04) who used the FORS instrument on the 
VLT to obtain radial velocities for 18 NGC 3379 inner GCs. 
Ten of the clusters in the P04 dataset are in our sample, in- 
cluding nine of our class A objects and one class B cluster. 
The agreement between the two studies is excellent with an 
average difference of 1 kms ' and a dispersion of 30 kms 
Since 30 kms ' is close to the typical uncertainty in the lower 
spectral resolution P04 data, this comparison is consistent with 
the small uncertainties we find for our velocities of only about 
10 km s"' in most cases. 



2.3. Separation of globular clusters and stars 

The histogram of measured radial velocities for the objects 
we observed is shown in Fig.|2 The heliocentric velocities of 
NGC 3379 and NGC 3384 are 904 km s ' and 721 kms ' re- 



http : //atlas . obs-hp . £r/elodie 



spectively (Falco et al. 1 1 9991 Smith et aI. I2000l l and the galax- 
ies are expected to have modest velocity dispersions. As a re- 
sult, GCs are well separated from stars in the data: there is a 
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200 km s"' gap between the lowest velocity object we consider 
a globular cluster and the highest velocity star. Specifically, 
we classify objects with 500 < u < 1300 kms"' as GCs, 
and objects with -102 < y < 332 kms"' as stars. The lat- 
ter range is in good agreement with the expected star field 
velocity of about 60 kms ' in the direction of NGC 3379 at 
(l,b) = (233?5,57?6). 

We then have a total sample of 42 GCs with radial veloci- 
ties, listed in TableQlalong with various properties; 34 of these 
GCs have highly confident class A radial velocities, and 8 have 
likely class B determinations. These form the basis for the anal- 
ysis and discussion in subsequent sections. 

We also have radial velocities for 36 stars in our field. This 
is consistent with expectations from the photometric survey of 
RZ04, which predicts that about 60% of our objects at a dis- 
tance of 4' from NGC 3379 are GCs, in agreement with the 22 
GCs and eleven non-GCs (stars and unclassified objects) within 
this radius. Most unclassified targets have V ^ 21 and are prob- 
able galaxies, compact but unresolved in the MOSAIC images 
(typical seeing r.'3) - including [On] emission line objects at 
0.35 0.55 for the LR4 wavelength range. The total num- 
ber of stars and background objects found is consistent with 
the work of RZ04, given our fibre allocation and V<22 limit. 

3. Properties of the globular cluster systems of 
NGC 3379 and NGC 3384 

3.1. Identifying globular clusters with NGC 3379 and 
NGC 3384 

The core of the Leo I group has a compact appearance on the 
sky with the central elliptical (El) galaxy NGC 3379 sepa- 
rated by only 7' from the nearby SBO galaxy NGC 3384, and 
by 10' to the small spiral NGC 3389, corresponding to pro- 
jected distances of 20 and 29 kpc respectively. However, the 
galaxies are fairly clearly delineated in their radial velocities, 
with NGC 3379 at 904 kms"', NGC 3384 at 721 kms"', and 
NGC 3389 with 1300 kms '. The velocity for NGC 3379 is 
the average of the Updated Zwicky Catalogue (UZC, Falco et 
al. ll999> and Smith et al. J2000> values, and the velocities for 
NGC 3384 and NGC 3389 are from the UZC. The scatter of 
measurements in the literature suggests uncertainties of 10- 
20 km s"' in these numbers. 

In Fig. |3lwe show the location and velocities of our GC 
sample overlayed on an image of the Leo I group. This fig- 
ure shows the clear velocity offset of the group of GCs around 
NGC 3379, which all scatter around the galaxy's systemic ve- 
locity, from the group around NGC 3384, which have veloci- 
ties like that of NGC 3384 and are thus offset to lower v than 
NGC 3379 and its GC system. The velocity difference of these 
two galaxies and their GC systems is useful for determining 
which clusters in their outer halos belong to which galaxies 
given their closeness in projection on the sky. 

Based on their radial velocities and positions, we assign 
all of the GCs spatially centered around NGC 3379 to its GC 
system, along with the GCs located in intermediate spatial re- 
gions that have velocities greater than the systemic velocity 
of NGC 3379, and the globular gc756. We also assign the 



eight GCs centered on NGC 3384 to this galaxy, and note 
that these have a strong rotation signature, which we discuss 
further in Sect. 13.31 We consider the objects gc573, gc571 , 
ad1102, and ad 1021 as possibly belonging to either galaxy. 
We note that it would seem unlikely that NGC 3379 has only 
outer globulars to its south and not its north, so some of these 
four would seem likely to belong to NGC 3379. On the other 
hand they are slightly closer in projection to NGC 3384, so 
some of the four would appear likely to be associated with its 
halo. The velocities of these four GCs are also generally con- 
sistent with NGC 3384 (see Fig.|4}, but if none of them is as- 
signed to NGC 3379, then its outer halo would have a strange 
bias towards only globulars with velocities greater than both 
the systemic velocity of NGC 3379 and its inner GC system 
(a similar argument supports the identification of gc756 with 
NGC 3379). We also note that although at very large radii it 
may become more sensible to consider the halo of the group of 
galaxies together, given the clear velocity offset of the galaxies 
and their inner GC systems, a closer relationship to the individ- 
ual galaxies is indicated over our field of view of x 30 kpc. 

For the calculations that follow, we use the assignments de- 
scribed above and given in Tabled and also consider the effect 
on these calculations of including or excluding various combi- 
nations of the GCs with uncertain assignments. As noted below, 
we find that the conclusions are not significantly affected by the 
choices in these assignments. 

3.2. Kinematics of the globular cluster system of 
NGC 3379 

Our sample of NGC 3379 GCs with radial velocities consists of 
the 30 clusters clearly identified with NGC 3379 and the four 
for which this assignment is uncertain. To these, we add eight 
additional GCs from the study of P04 (see Sect.l2.2>. 

For the sample of 38 GCs associated with NGC 3379, we 
find a mean of y = 947 km s ' with an uncertainty of 38 km s ' 
determined via Monte Carlo simulations. The projected veloc- 
ity dispersion of the full sample is o-p = 169;^,^ kms"'. We 
discuss the o-p as a function of radius in the following sec- 
tion, in which we use this to constrain the mass of the halo 
of NGC 3379. One additional note is that the mean value of the 
GC radial velocities is Icr larger than the velocity of the galaxy 
at 904 kms ' (Sect. [Till- Our GIRAFFE data have good spec- 
tral resolution and are highly repeatable from night to night. 
Moreover, there is excellent agreement between our radial ve- 
locities and the independent data of P04. Thus, the most likely 
explanation is that the small offset is a statistical fluctuation, 
which can be tested with more data. For modelling purposes, 
we will use hereafter the mean GC velocity of 937 km s"' from 
the central regions with good azimuthal coverage (Fig.l^. 

Our data also show little evidence for rotation in the 
NGC 3379 GC system. Although our sample size is modest, 
even with a much smaller sample we are able to detect rota- 
tion in NGC 3384 owing to our small velocity uncertainties 
(see Sect, \3.3l . The exact constraint on the rotation around 
NGC 3379 depends somewhat on the sample definition, but by 
comparison to the strong signal seen around NGC 3384, the 
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Table 1. FLAMES heliocentic radial velocities v of globular clusters in Leo I: for each of the 42 confirmed GCs, v and its error 
estimated by fxcor are the median values over the ELODIE templates (ordered by projected distance to NGC 3379, R in arcsec; 
the position angle is counted in degrees from to 90, N to E). Numbering is from Rhode & Zepf J2004> as well as the V 
magnitudes, B -V and V - R colours. Equatorial J2000.0 coordinates should be accurate to < 0'.'2 with respect to the GSC2.2. 
The Tonry & Davis ( I1979> cross-correlation coefficient "R^^^ is the median value over the ELODIE templates. Class B objects are 
commented, as well as the 10 GCs in common with Puzia et al. (I2004> with their velocity estimate. Full tables of all MEDUSA 
targets - including stars and galaxies/unclassified objects - are available online. 
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" = Brodie & Larsen jMSl N3384-FF-7 "fuzzy" cluster. They found u = 905 ± 23 km s"' (red channel value) using A:ecA:/LRIS. 



overall rotation in the NGC 3379 GC system appears to be less 
than 100 km s ' about any axis. This low rotation is consistent 
with that seen in the stars at small radii (e.g. Statler & Smecker- 
Hane 1 199*^1 and is similar to that seen in the GC system of the 
elfiptical galaxy NGC 4472 (see Zepf et al. l2000l Cote et al. 
pOHT) although probably less than that in NGC 5128 (Peng et 
al. .2004j . a peculiar merger remnant. 



3.3. Kinematics of tlie globular cluster system of 
NGC 3384 

Our sample also has eight GCs that are clearly associated with 
NGC 3384. Even with just these eight radial velocities, there is 
evidence for rotation {v/cr ~ 2) in the GC system of NGC 3384, 
which can be seen visually in Fig.|3 Specifically, the four GCs 
to the south-west of the galaxy centre have radial velocities 
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Fig. 3. Velocities of the 42 confirmed GCs overplotted on a V- 
band MOSAIC sub-image of 25' x 15' (north is up, east to the 
left). Crosses (x) are GCs with v > 904 km s ' which we adopt 
as the systemic velocity of NGC 3379. Circles (O) are those 
GCs with V < 904 kms the symbol size is proportional to 
the offset, up to about ±400 kms"', with class A objects in 
boldface. Note the clear rotation of the eight GCs along the 
NGC 3384 disc and their lower typical velocities compatible 
with the systemic value of this SBO. The seven outer (R > 5') 
GCs associated with NGC 3379 and the four globulars with 
uncertain assignments to one of the two galaxies are identified 
individually, with the unassigned objects in italics. 



200 kms"' lower than the four to the north-east. This rota- 
tion signature is very similar to that found by a study of 68 
PNe around this galaxy by Tremblay et al. J1995> . The rota- 
tion in the GC system of NGC 3384 suggests that rotation may 
not be unusual in the GC systems of disc galaxies. This brings 
into question whether the lack of rotation in the large majority 
of the Milky Way globular cluster system (e.g. ZinnQ985j is 
somewhat unusual, and the rotation seen in the Messier 3 1 GC 
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system (e.g., Perrett et al. l2002t is possibly more typical. 



Fig. 4. Relative velocities vs. distances to NGC 3379 for the 
spectroscopic sample of 42 globular clusters. GCs associated 
with NGC 3379 are shown as squares (bold for class A; di- 
amonds are eight additional GCs from P04), those associated 
with NGC 3384 as filled circles and GCs possibly associated 
with either galaxy as open circles. 

3.4. Colours of the globular clusters around NGC 3379 
arid NGC 3384 

We can use our dataset to examine questions about the colours 
of the NGC 3379 and NGC 3384 GCs with a specti-oscopically 
confirmed dataset. In Fig.|5lwe plot the B — R colours from 
RZ04 of our GC sample as a function of their distance from the 
centre of NGC 3379. The data show that NGC 3379 has more 
blue, metal-poor GCs than red metal-rich ones, in agreement 
with the suggestion of RZ04 from their photometry. As noted 
in that paper, of current models for the formation of GC sys- 
tems, only those that involve galaxy mergers (e.g. Ashman & 
Zepf'1992J Beasley et al. '2002i are readily able to account for 
this larger number of blue relative to red GCs in an elliptical 
galaxy like NGC 3379 with its modest luminosity and small 
GC population. 

The data shown in Fig. |5] also indicate a colour gradient 
in the GC system of NGC 3379, such that the ratio of blue to 
red GCs increases with distance away from the galaxy centre. 
This is consistent with RZ04 and many other previous stud- 
ies for the globular cluster systems of a number of elliptical 
galaxies (see Ashman & Zepf 1199"^ for a review). However, 
RZ04 also noted that the photometry indicated that there may 
be red, metal-rich GCs at large distances around some ellip- 
tical galaxies. In NGC 3379, our data show a metal-rich GC 
(gc581) at a projected distance of almost 6', corresponding to 
about 17 kpc. The GC ad1 021 is also interesting in this regard, 
as it is very red, and is almost 10' away from NGC 3379 and 
about 5' away from NGC 3384. Thus there are two spectro- 
scopically confirmed, metal-rich clusters well into the halos of 
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Fig. 5. B-R colours vs. radii for our spectroscopic sample of 
42 GCs (symbols are as in Fig.l^. The plot shows the trend of 
redder colours for GCs close to the NGC 3379 centre and bluer 
colours at large radii. 

these galaxies. We also note that two of the inner GCs associ- 
ated with NGC 3384 are very red, which might be attributable 
to dust in the system. 




R (arcmin) 

Fig. 6. Projected velocity dispersion vs. radius for the 
NGC 3379 GCs. The error bars show radially-binned data, 
where the vertical bars are for the +34% uncertainties, and the 
horizontal bars show the range covered by 68% of the data. 
The dotted line shows a moving window estimate, with dashed 
lines signifying the ±34% uncertainties. The solid line shows 
the best-fit scale-free power-law profile. The boxes at the bot- 
tom indicate the radial positions of the GC velocity data. 



4. Constraints on the halo of NGC 3379 

4.1. Velocity dispersion profile 

One of the primary goals of this study is to determine the radial 
velocities of GCs to probe the halo dynamics of NGC 3379. To 
test mass models of the halo of NGC 3379, we use our dataset 
to determine the projected velocity dispersion at several dif- 
ferent radii. As shown in Fig.|6j we adopt three radial bins to 
allow the comparison of the mass over a range of radii, corre- 
sponding in the inner region to the largest radial extent of the 
stellar data, in the intermediate region to the largest extent of 
the PNe data, and the outer region which is uniquely probed 
by our globular cluster dataset. We then plot the projected ve- 
locity dispersions o-p in these three radial regions in Fig.|6] To 
calculate o-p, the mean velocity of uq - 937 kms"' is adopted, 
and the measurement uncertainties 6ui are taken into account. 
The +34% uncertainties are estimated by 30 000 Monte Carlo 
simulated data sets drawn from the best-fit model. Note that the 
dispersion o-p discussed throughout this paper incorporates any 
rotational support, and thus is technically Ui-^^ = (fp + o"p)'^~- 

Figure |6l shows the binned velocity dispersion profile. We 
estimate o-p = 205+^^ kms"', 155!^^ kms"' and 147+^4 kms"' 
at R~0'.8 (13 GCs), -2:5 (18 GCs) and -11' (7 GCs), respec- 
tively. Inclusion of the four uncertain GCs has very little effect: 
the outer o-p becomes 147;^ kms"^ Alternatively adopting 
904 km s"' for the NGC 3379 velocity has the effect of raising 
the outer two o-p estimates to 165;|;27 ^ ' l^S;!;^^ km s"' . 



We also construct a smoothed estimate of (Tp{R), using a 
"moving window". This performs a o-p estimate at every radius 
R using the nearby data (each estimate is thus not independent). 
Unlike previous versions of this method which used a window 
of constant width (Zepf et al. 120001 Cote et al. 1200 U or with 
fixed number of data points (Kissler-Patig & Gebhardt ll9981 
Cohen EDOOI I. we use an exponential weighting function for the 
contribution of each data point v, to the estimate at radius R. 



Wj{R) = — exp 
o-R 



iR-Rf 



(1) 



where the half-width ctr oc [/?2(/?)] and l.(R) is the underly- 
ing surface density distribution from which the measurements 
are drawn. For "perfect" sampling, the number of points con- 
tributing to each estimate is constant. For optimal smoothing, 
the number of points within +o-« has a maximum of ~ 25. For 
the current data set, the surface density is the R^^'^ law found by 
RZ04, and the window width is normalised to o-^ ^ 1 ^8. 

Figure|6]shows the result of this procedure. The dispersion 
profile is consistent with being constant over the wide range 
of radii studied, with a possible increase to small radii. The 
inclusion of the four uncertain GCs does not affect the pro- 
file much, but a mean velocity of 904 kms"' would raise o-p 
in most places by 10-20 kms"'. Utilising a maximum likeli- 
hood technique for fitting the data to a scale-free power-law 
(see Appendix A), we also find the same result that the veloc- 



8 



G. Bergond et al.: Wide-field kinematics of globular clusters in the Leo I group 



ity dispersion is generally constant with radius with a possible 
modest increase in the dispersion to small radii. 

4.2. Halo models 

A primary goal of this programme is to test whether our 
measured velocity dispersion profile for the GC system of 
NGC 3379 is consistent with the mass distribution expected 
from standard ACDM halos or other models. To perform this 
test, we construct mass models of NGC 3379 with and without 
dark matter halos. The circular velocity of these models is: 



+ y: 



r \ rj r + r. 



(2) 



where and a* describe the stellar mass distribution (after 
Hernquist I 19901 1 and and r, are the Navarro, Frenk & White 
J1996> NFW parameters of the dark matter halo. More accurate 
representations of the stellar and dark matter distributions are 



possible (e.g. Mamon & Lokas 120051 . but at this stage we wish 
to use identical models to R03 for consistent comparisons. 

For the stellar mass distribution, extensive studies of stel- 
lar dynamics in the inner regions of NGC 3379 (Kronawitter 
et al. 



120001 Gebhardt et al. 120001 Cappellari et al.|^05|, find 
stellar mass-to-light ratios of Y*3 - 5-6 (in units of Tq b). 
These models neglect any dark matter found inside Rg, so the 
actual T*,B may be lower The orbit models of R03 separate 
M{r) into the stellar and dark components, giving T* ^ = 6 
(based primarily on the kinematics of the stars rather than the 
PNe). Stellar population synthesis models with a Kroupa IMF 
imply Y^,,B = 4-7 (Gerhard et al. 2001:. Napolitano et al. 200j 
Cappellari et al. I2005I I. We adopt ^ = 5.8 as the overall 
best estimate, and combine this with the observed light profile 
to construct the galaxy mass profile. The resulting stellar mass 
models for NGC 3379 gives = 603 kms ' and a* = 03 
(^e - 0'6) for use in Eq. (|2j above. 

To determine the and values for a given ACDM halo, 
we compute the virial mass and radius, My-„ and ryii , and ex- 
press the model halos in terms of Myii and the concentration 
Cvir = rvii/rs. In ACDM, there is a well-determined mean rela- 
tion between these parameters: Cvir - 17.9(Mvir/10"MQ)"'''^^, 



with a 68% log scatter of 0. 14 (Bullock et al. l200H Napohtano 



et al. l2005t . Our approach here differs from R03 in assuming a 
fluctuation amplitude erg = 0.9 rather than 1 .0, and in using a 
virial overdensity of 101 rather than 200 - changes that permit 
less massive inner halos in ACDM. 

We then consider whether ACDM halos with various 
masses and concentrations are consistent with the velocity dis- 
persions we find for the NGC 3379 GCs. To do this, we cal- 
culate the projected velocity dispersion of the globular cluster 
system using a Jeans equation: 



2 r°° rdr r 



j{r')vl{r')dr' 



r' 



(3) 



In making this calculation of the crp(R) expected for the GC 
system given a mass model for the galaxy, we assume the fol- 
lowing: 1) spherical symmetry; 2) isotropic distribution func- 
tion and Gaussian LOSVDs; 3) density distribution of tracers 



j(r) well-known; 4) dynamically negligible gas content; 5) stel- 
lar mass-to-light ratio Ti, well-constrained by other techniques. 

The spherical assumption is motivated by the roundness 
of NGC 3379 (e.g. Peletier et al. .1990J , the absence of evi- 
dence for substantial flattening of the GC system of NGC 3379 
(RZ04), and the evidence that globular cluster systems appear 
to be typically at least as round as their host galaxies (e.g.. 



Ashman & Zepf lI998l l. Similarly, isotropy is observed for sys- 



tems where sufficient data are available (e.g. Zepf et al. 120001 
Cote et al. i2003t ; our a-p fits to the NGC 3379 GC velocities 
give very similar results to traditional cTp calculations, indicat- 
ing that Gaussianity is a good approximation. The density dis- 
tribution we adopt is a Hernquist (.1990 J model which matches 
the jjj^ flttejj to the NGC 3379 globular clusters (RZ04): 
j(r) oc r\r + fl)-^ where a = 9^8 {R^ = 17;9). The mass of 
hot gas within 100 kpc of NGC 3379 is no more than ^ IO'^Mq 
(David et al. 120051 conservatively assuming a linear increase in 
gas mass from the inner parts), which is a negligible fraction of 
the stellar mass of ^ IO^Mq. 



4.3. Comparison of globular cluster system velocity 
dispersion to models 

The comparison of the o-p{R) of the GC data to predictions for 
various models of the mass distribution for NGC 3379 is shown 
in Fig.0 The GC data are in good agreement with a wide va- 
riety of ACDM halos that have typical concentrations for their 
masses and match the inner stellar data. In detail, the range of 
ACDM halos plotted is based on halos which fall within the Icr 
scatter in the halo mass-concentration relation of Bullock et al. 
J2001I I, with the additional constraints that the stellar mass to 
dark halo mass is not lower than typical for galaxy clusters and 
is not larger than the cosmic baryon fraction. 

A valuable constraint on the halo mass of NGC 3379 is 
provided by observations of a rotating ring of H i gas at R ~ 
100 kpc (Schneider et al. [T^ Schneider [1985}. The posi- 
tions and velocities of the Hi are consistent with the gas be- 
ing in a Keplerian orbit around the barycentre of NGC 3379 
and NGC 3384, with about of the mass within 100 kpc be- 
longing to NGC 3379. This then gives a mass for NGC 3379 
of about 4-5 xIO^Mq at R ~ 100 kpc and Yb - 25-35. 
A ACDM halo model that matches the Hi ring and is gener- 
ally consistent with the GC and FN dynamical constraints is 
shown in Fig.0 This "consensus" model predicts a GC veloc- 
ity dispersion that is close to that observed, with a preference 
for smaller values. 

This "best-fitting" ACDM halo is only slightly below the 
GC data and has a concentration well within the normal range 
for a halo of its mass. The only notable feature of this ACDM 
halo is that it gives a fairly small mass-to-light ratio, with 
Tb - 35. Modest mass-to-light ratios for early-type galaxies 
with luminosities around L* like NGC 3379 may be in agree- 
ment with expectations from some recent halo occupation mod- 
els (e.g. van den Bosch et al. I2003> . although the value for 
NGC 3379 is about a factor of two less than the model aver- 
age for galaxies of its mass. Another implication of the inferred 
mass-to-light ratio is that if the halo has the universal baryonic 
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fraction of 0.17 of the concordance model, and given the 
of about 6 measured in the inner regions, the for the halo 
would suggest most of the baryons have been turned into stars 
(see Napolitano et al.;2005 for this calculation for a number of 
galaxies). 

Our data also allow us to assess alternatives to ACDM ha- 
los. Assuming no dark matter and a Newtonian force law, we 
expect o-p = 98 kms ' at R ^ I'.S and 49 kms ' at 10'. This 
mass traces light model is clearly ruled out by the data (see 
Fig.Q. We also consider MOND, which has been proposed to 
provide a good fit to the earlier PNe data (Milgrom & Sanders 
I2003> . The comparison of our observed o-p with the predic- 
tions of MOND is shown in Fig.0 The velocity dispersion of 
the GCs is slightly higher than the MOND prediction, although 
the uncertainties mitigate against an unambiguous rejection of 
MOND. 

A final comparison is to the published velocity dispersions 
of NGC 3379 PNe (R03). The reasonable agreement of all of 
the constraints with the overall "best-fit" mass model indicates 
that the GCs and PNe are broadly consistent. In detail, although 
the PNe have a cTp of about 120 kms"' in the radial range of 
on < R < 4' where the PN data overlap with our GC data, 
compared to a o-p of about 170 km s ' for the GCs in the same 
region, the two do not give greatly discrepant masses. This is 
partly due to radial anisotropy in the PNe system suggested by 
orbital modeling of the PN system (R03), which naturally gives 
a somewhat lower o-p for the same mass, and partly due to the 
uncertainties in o-p from the modest number of radial velocities. 
Furthermore, the GCs appear to have a shallower spatial profile 
than the galaxy light (and by inference the PNe) in the relevant 
region (RZ04) which would produce a smaller o-p for the PNe 
relative to the GCs. If additional data find a difference in the o-p 
beyond these effects, several possibilities present themselves. 
These include the possibility that the GCs may have modestly 
tangential orbits, that orbits may be disturbed by group inter- 
actions, or that NGC 3379 may have a substantial disc compo- 
nent seen face-on (e.g. Capaccioli et aI. II991> which would be 
more likely to include PNe than GCs. The question of the halo 
mass profile inferred for NGC 3379 has also been addressed 
in a recent paper by Dekel et al. J2005I I. They simulated disk 
galaxy mergers and found that in the resulting elliptical galax- 
ies, tracers like PNe have steeply declining o-p profiles because 
of radially biased orbits, flattened triaxial structures, and steep 
spatial profiles. They predicted that GCs would be less prone to 
these effects and should have flatter dispersion profiles, which 
is consistent with our results. 

5. Conclusions 

We have used the FLAMES multi-fibre spectrograph on the 
VLT to obtain radial velocities of globular clusters (GCs) 
around the Leo I group galaxies NGC 3379 and NGC 3384. 
FLAMES allowed us to obtain spectra of 42 GCs over a wide 
field (V ^ r ^ 14') with very small uncertainties (~ 10 kms"'). 

1 . Our primary conclusion is that the projected velocity dis- 
persion of the GC system of NGC 3379 is consistent with that 
expected from standard dark matter halos in ACDM cosmolo- 
gies. This is based on comparing our observed velocity disper- 
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Fig. 7. Model predictions compared to the observed velocity 
dispersion profile of the NGC 3379 GCs. Error bars are as de- 
scribed in Fig.|6l The hatched region describes profiles for typ- 
ical ACDM halos as described in the text. The solid line shows 
the "consensus" model that matches the H i data and provides 
a reasonable fit to the GC and PNe results. The dotted line rep- 
resents a mass traces light model with Newtonian dynamics, 
while the dot-dashed line is with MOND. 

sions for the GC system over a wide range of radii to the expec- 
tation from the Jeans equation for the mass profile of ACDM 
halo models and the observed radial profile of the GC system. 

2. We find evidence that the GC system of the SBO galaxy 
NGC 3384 has a rotation similar to that of the stellar disc and 
planetary nebula system of this galaxy. This result suggests 
that significant rotation in the globular cluster systems of disc 
galaxies may not be rare. 

3. We find a colour gradient in our spectroscopic sample 
such that the fraction of red, metal-rich GCs decreases from 
small to large radii, in agreement with photometric studies. 
Even with this colour gradient, we also find some red GCs at 
large distances from the centre of their host galaxy, which may 
place interesting constraints on models for the origin of metal- 
rich globular clusters. 
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Appendix A: Fitting dispersion profiles to discrete 
velocity data 

Radial binning or smoothing is not ideal for characterising data 
and testing models, as the grouping of data discards useful in- 
formation about the radial locations Rj and may skew the result 
by assuming cTp is constant within each bin. As an alternative, 
we consider a maximum-likelihood technique for fitting dis- 
crete velocity data to local velocity dispersions (see Ciardullo 
et al. 119931 Saglia et al. EOOOj . This technique allows one to fit 
data to a model for crp{R), which may be e.g. a constant disper- 
sion, a power-law, or the prediction of a Jeans model. 

We convolve a Gaussian line-of-sight velocity distribution 
(LOSVD) with a Gaussian measurement function y,- ± Svj to 
derive the likelihood of a given measurement y,- given a model 
dispersion crp(i?,); 



1 



: exp < , , 

which can be converted into ax^ statistic: 



(A.l) 



X 



Zj (vi-vo)- , r 2 x^i 



(A.2) 



This method assumes Gaussian LOSVDs, so strictly speaking, 
we are computing the best-fit Gaussian parameter tXp rather 
than the true second-moment o-p. For near-isotropic systems, 

O-p - CTp. 

To apply this to our globular cluster data for NGC 3379, 
we characterise the dispersion data by a scale-free power-law 
profile: 



(A.3) 



where Rq = 2'.2. We find ctq = 168!^i kms"' and = 0.25 + 
0.25 (see Fig. |6}. This result is scarcely affected by including 
the four uncertain clusters; on the other hand, adopting a mean 
velocity of 904 kms"' would mean ctq - lll^J^^ kms"' and 
rp = 0.20 + 0.25. 

This dispersion estimator is checked using Monte Carlo 
simulations. A small correction for estimator bias is derived 
using the difference between the simulations' input and median 
output. The (doubly-debiassed) 68% scatter of the simulations 
gives the likely range of the dispersion. 
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Table A.l. Heliocentric radial velocities v (in km s ' ) of FLAMES candidate globular clusters in Leo I. The first 34 ("Class A", 
ordered by V-band magnitude), are certain; the next 8 ("Class B") are probable; the final 8 ("Class C") are possible. Object ID 
and photometry are from Rhode & Zepf J2004> - the "ad" prefixed objects have been added by hand to their initial list of GC 
candidates, h^^^^, vv^-cf and H^^^ are the height, width (in pixels) and Tonry & Davis ( 1979: ) quality coefficient of the Gaussian fit 
to the cross-correlation function peak. The inferred fxcor velocities and their estimated errors are the median values over the 
ELODIE templates. 
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Table A.2. Heliocentric radial velocities v (in km s ' ) of FLAMES candidate stars. The first 32 ("Class A"), ordered by V-band 
magnitude, are certain; the next 4 ("Class B") are probable; the final one ("Class C") is a possible star Object ID and photometry 
from Rhode & Zepf J2004> . h^^^, vv^^p and "R^^p are the height, width (in pixels) and Tonry & Davis J1979> quality coeflicient 
of the Gaussian fit to the cross-correlation function peak. The inferred fxcor radial velocities and their estimated errors are the 
median values over the ELODIE templates. 
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Table A.3. The 21 unclassified FLAMES targets, ordered by V-batid magnitude. Object ID and photometry are from Rhode & 
Zepf (|in03- 
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gel 080 


10:47:30.38 


+ 12:40:53.4 


21.16 


0.73 


0.37 




gc970 


10:47:38.84 


+ 12:44:27.6 


21.19 


0.53 


0.39 




gc614 


10:48:03.92 


+ 12:39:46.4 


21.43 


1.06 


0.48 




gc873 


10:47:45.16 


+ 12:33:17.7 


21.62 


0.66 


0.46 




ad1106 


10:48:10.08 


+ 12:24:08.9 


21.64 


0.73 


0.40 


[On] at; -0.40 


gc718 


10:47:53.88 


+ 12:43:16.6 


21.66 


0.52 


0.42 




gc511 


10:48:12.27 


+ 12:27:42.2 


21.70 


0.76 


0.38 




gc803 


10:47:48.84 


+ 12:45:03.2 


21.74 


1.04 


0.64 




gel 075 


10:47:30.95 


+ 12:31:06.6 


21.74 


0.84 


0.49 




gc738 


10:47:52.79 


+ 12:27:27.7 


21.76 


0.91 


0.50 


Emission lines? 


gc849 


10:47:46.57 


+ 12:42:35.9 


21.78 


0.82 


0.47 


[On] at z~ 0.38 


gc983 


10:47:37.91 


+ 12:35:11.8 


21.78 


0.82 


0.50 




ad697 


10:48:25.26 


+ 12:39:11.8 


21.82 


0.94 


0.50 




gc292 


10:48:28.21 


+ 12:38:28.6 


21.84 


0.64 


0.42 




gc419 


10:48:18.12 


+ 12:34:41.9 


21.84 


0.90 


0.45 




gc344 


10:48:23.74 


+ 12:29:17.9 


21.93 


0.69 


0.45 


[On] at z~ 0.40 


gc600 


10:48:04.85 


+ 12:35:38.3 


21.93 


0.77 


0.48 





